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ABSTRACT 

The 21cm-galaxy cross-power spectrum is expected to be one of the promising probes 
of the Epoch of Reionization (EoR), as it could offer information about the progress 
of reionization and the typical scale of ionized regions at different redshifts. With 
upcoming observations of 21cm emission from the EoR with the Low Frequency Array 
(LOFAR), and of high redshift Lya emitters (LAEs) with Subaru’s Hyper Suprime 
Cam (HSC), we investigate the observability of such cross-power spectrum with these 
two instruments, which are both planning to observe the ELAIS-Nl field at z = 
6.6. In this paper we use N-body + radiative transfer (both for continuum and Lya 
photons) simulations at redshift 6.68, 7.06 and 7.3 to compute the 3D theoretical 21cm- 
galaxy cross-power spectrum, as well as to predict the 2D 21cm-galaxy cross-power 
spectrum expected to be observed by LOFAR and HSC. Once noise and projection 
effects are accounted for, our predictions of the 21cm-galaxy cross-power spectrum 
show clear anti-correlation on scales larger than ^ 60h~^ Mpc (corresponding to 
k ^ O.l/i Mpc“^), with levels of significance p = 0.04 at z = 6.6 and p = 0.048 at 
z = 7.3. On smaller scales, instead, the signal is completely contaminated. 

Key words: galaxies: high redshift - cosmology:observations - reionization - inter- 
galactic medium 


1 INTRODUCTION 

The Epoch of Reionization (EoR) is one of the greatest ob¬ 
servational frontiers in modern astrophysics. It corresponds 
to the transition from a neutral to an ionized Universe, as 
mostly young, star-forming galaxies reionized the intergalac- 
tic medium (IGM) surrounding them. 

The early stages of cosmic reionization are character¬ 
ized by the presence of isolated HII regions, which grow in 
number and size as structure formation progresses. During 
the latest stages of the process, the ionized regions are very 
large, with individual sizes exceeding tens of comoving Mpc. 
Eventually, they overlap and IGM reionization is complete 
(e.g. Gnedin fc Ostriker| 1997[ Giardi et al.||2000[ Gnedin 


|2000|[Ciardi, Ferrara fc White|2003||Maselli, Ferrara fc Ci^ 


rdi|20031 Santos||2004| Furlanetto, Zaldarriaga & Hernquist 

2004 

Haiman & Gen|2005||Iliev et al.|20061 McQuinn et al. 

2007 

Zahn et al.|2011 Giardi et al.|2012 Iliev et al.|2014|. 

While early simulations were only a few Mpc in size ( 

Gnedin 

& Ostriker 1997 Giardi et al. 2000 Gnedin 2000 

1, novel 


codes for cosmological N-body and hydrodynamical simu¬ 
lations and for radiative transfer finally enable reionization 
simulations with volumes larger than ~ 100 Mpc ( Iliev et al. 
2006 McQuinn et al.|2007 Iliev et ^|2014 1, allowing for a 


correct abundance of rare massive haloes (Barkana & Loeb 


2004 Li et al. 2007 Trac & Gnedin 20111, while resolv¬ 
ing also dwarf-size galaxies with masses ~ 10® M©, which 
are considered to be the main sources of ionizing photons 
(Volonteri fc Gnedin|2009 Loeb|2009 Robertson et al. 2010 


E-mail: dvrbanec@mpa-garching.mpg.de 


Fontanot, Gristiani fc Vanzella|2012 1. In recent years, a sig¬ 
nificant improvement has taken place in terms of dynamic 
range of such simulations and sophistication in the imple- 
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mentation of various relevant physical processes (e.g. Trac 
fc Gnedi^|2011 1 . 


Absorption spectra of high-redshift quasars suggest that 
reionization was completed by 2 « 6 ([Fan et al.||2006[ 


Bolton et al.||20ll McGreer, Mesinger & D’Odorico 20151. 

On the other hand, measurements of the primordial Cos¬ 
mic Microwave Background (CMB) radiation obtained by 
the WMAP satellite indicate that the process started much 
earlier, suggesting that the Universe was neutral until 2 = 


lO.lil.O, if instantaneous reionization is assumed (Komatsu 
et al.| 2011 |. More recent measurements with Planck instead 


give 2 = 8.8±1.1 ( Planck Collaboration et al.|2015 1. Present 
observations though do not offer much information neither 
on the progress of reionization nor on the main sources re¬ 
sponsible for it. Detection of the 21cm line from neutral 
hydrogen promises to offer insight in this respect. 

Because of the wealth of Information such observations 
could provide, there are significant efforts to detect reion¬ 
ization by mapping the 21 cm hyperfine transition lin e of 
neutral hydrogen w ith radio arrays such as LOFAl-0 (|vm| 


Haarlem et al. 


2013| ), MWA[5 PAPEE[5 CMRlj^and SKI 


Calculations predict that the cosmological 21cm signal from 
the EoR will be extremely faint, while the system noise 
and the foregrounds will be orders of magnitude larger (e.g. 


Madau, Meiksin & Rees(1997| 

Ciardi & Madau|2003t Zaldar- 

riaga, Furlanetto & Hernqulst 

20041 1’urlanetto, Oh & Briggs 

20061 McQuinn et al. [20061 Mellema et al. [12006 IMorales, 


Bernard! et al.|2009[ Pober et al.|2013| Parsons et al. 12014 1. 

Due to the low signal to noise ratio (which for LOFAR is 


Labropoulos et al. [2009j^ Harkere^^aD 201^ G«il^^^Gaensl^ 

fc Wyithe||2011[ Chapman et al.|2012| ‘[^aroiibi et al.||2012[ 


0.2, Labropoulos et al. 2009 Zaroubi et al. 2012 1 the 


Hrst observations of the 21 cm signal will measure only sta¬ 
tistical properties, such as the rms and power spectrum of 
the brightness temperature and their evolution with time 
Ciardi fc Madau] 2003| Morales &: Hewitt] 2004[ Barkana 


Loeb| 2005jJBowmaipMorales^^^HewittJ|2006|jMcQuim 

et al. 2(306| ^itchard fc FurlanettopOOfl' Jellc et al.||2008[ 


Pritchard & Loeb 2008| Harker et al. 2009[ 20101. Other 
probes of reionization that will be possible with observa¬ 


tions of the redshifted 21cm line are imaging (e.g. Koop 
mans et al.||2015l [Mellema et al.|2015| [Wyithe, Ceil fc Kim 


20151, which will not be possible before SKA, 21cm forest 
] Mack fc Wyithe|2012 Ciardi et al.||201^ 20151, and cross¬ 
correlation with other observations at different wavelengths. 


yaev 


such as near-infrared background radiation (NIRB; e.g. Fer¬ 
nandez et al. 2014 Mao 20141, kinetic Sunyaev-Zel’dovich 
effect (kSZ; e.g. Jelib et al.|2010 |Tashiro et al.|2010[), galax- 
ies (e.g. Lidz et al. 2009 Wiersma et al. 2013 Park et al. 


20141, CO line (e.g. Righi, Hernandez-Monteagudo & Sun 
2008| Visbal fc Loeb|2010[ Lidz et al. 2011 1 and CH line 


Gong et al.||2012 Silva et al.||2015 1. Detection of such 


(e.g. _ 

cross-correlations can provide further insight into different 
aspects of the EoR, such as the progress of reionization and 
the redshift at which the process is halfway, the evolution 


^ http://www.lofar.org 

^ http://web.haystack.mit.edu/arrays/MWA 
^ http://eor.berkeley.edu 
^ http://gmrt.ncra.tifr.res.in 


of the neutral hydrogen content, and the typical scale of 
ionized regions at different redshifts. 

Another way to explore reionization is to probe high- 
2 , young, star-forming galaxies, which are considered to be 
the dominant sources of ionizing photons. Such galaxies are 
expected to have a strong Lya emission line due to the in¬ 
teraction of the interstellar medium with ionizing radiation 
from young massive stars ( Partridge fc Peebles|[T9^ I. De¬ 
pending on the detection method, such galaxies are typi¬ 
cally referred to as Lya emitters (LAEs) and Lyman break 
galaxies (LBGs). Star-forming galaxies that are luminous 
enough to be detected with existing telescopes most likely 
populate fairly mass ive dark matte r haloes, with masses in 
excess of 10^° Mq (Dijkstra 2014 b They ionize their sur¬ 


roundings forming large HII bubbles in which one or more 


star-forming galaxies reside (e.g. Dijkstra 2014 b Lya pho¬ 
tons emitted by those galaxies can therefore propagate and 
redshift away from line resonance through the ionized IGM 
before entering the neutral IGM | Miralda-Escude &: Ree^ 


]^998jJSantos|20d4[jCen^^&HairnaiJ2000^ Barton et al. 

Cen, Haiman fc Mesin ger||2005[ McQuinn et al. [2007j^ Iliev 

et al.||2008[ [Mesinger ST Furlanetto||2008| [Curtis-^Lake et al. 


2012 b These photons are then less likely to be scattered out 


of the line of sight. This is why LAE luminosity functions 


([Haiman fc Spaans|1999[ Malhotra & Rhoads|2004[ Haiman 

fc Cen[[20051 

Le Delliou et al.|2005 

IFurlanetto, Zaldarriaga 

fc Hernquistj 

20061 Le Delliou et a 

.[20061 Mao et al.[[2007 

Dijkstra, Wyithe fc Haiman|2007 

Jensen et al.|2013 20141, 

number density (Malhotra fc Rhoads[ 2006 1 and clustering 

(Furlanetto, 

Zaldarriaga & Hernquist|2006| McQuinn et al. 


methods to study the EoR with galaxies. A reduction in the 
number of observed sources, and thus a suppression of the 
luminosity function, is expected with increasing redshift, due 


to the larger amount of neutral gas in the IGM (e.g. Haiman 
fc Spaans|1999 b 

As of now, 207 LAEs have been observed at 2 = 6 . 45 - 
6.65 ( |Ouchi et al.|201^ , 1 at z=6.96 ( |Ota et al.|2008P and 7 
at 2 = 7.3 ( Konno et al.|[20l4 1, while there are no confirmed 
LAEs at higher redshifts. These numbers are likely to in¬ 
crease soon, as the Subaru telescope with its new Hyper 
Suprime-Cam (HSC) started narrow-band observations at 
redshift 6.6 and 7.3. The HSC has a 1.5 deg (in diameter) 
Field of View (FoV) and it will observe LAEs in its Deep 
( 2 = 6 . 6 , for a total of 28 deg^) and Ultra-deep ( 2 = 6.6 and 
7.3, for a total of 3.4 deg^ at each redshift) helds (Miyazaki 


et al. 2012 b It Is estimated that HSC will observe ~ 5500 
LAEs at 2 = 6.6 and ~ 40 LAEs at 2 = 7.3 (Ouchi 2012, pri¬ 
vate communication). 

Although detection of the 21cm signal and observations 
of LAEs will provide invaluable insight on reionization and 
its sources, the shape and normalization of their cross-power 


spectrum can offer ad ditional information ( Lidz et al.|2009 
Park et al. 2014[p] In this article we 


Wiersma et al. 


2013 


will present theoretical 3D 21cm-galaxy cross-power spectra 
computed from full radiative transfer + N body simulations 


by Iliev et al. (20141, post-processed with a Lya radiative 


transfer code by Jensen et al. (20131 and Laursen (2010 b 


http: / / www.skatelescope.org 


® [Wiersma et al.| ( |2013[ have shown that LAEs are more promis¬ 
ing probes of reionization than LBGs. 
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Differently from previous papers, here we accurately model 
LAEs, so that their simulated properties reproduce the ob¬ 
served ones. We will also show 2D 21cm-galaxy cross-power 
spectra from mock observations obtained by adding instru¬ 
mental effects for the LOFAR and Subaru telescopes. The 
paper is organized as follows: in section 2 we describe the 
simulations used; in section 3 we show the theoretical 3D 
cross-power spectra as well as the observed 2D cross-power 
spectra; in section 4 we show 21cm-galaxy cross-correlation 
functions. In section 5 we discuss the noise from both LO¬ 
FAR and HSC observations, and draw conclusions in section 
6 . 


The following set of cosmological parameters was used: 
Da = 0.73, = 0.27, Dt = 0.044, fe = 0.7 and ag = 


0.8, ris = 0.96, consistent with WMAP 5-year data (Ko- 
matsu et al.||2009f 


2 SIMULATIONS 


To compute the 21cm-galaxy cross-power spectrum we have 
used a full radiative transfer -f N-body simulation of reion¬ 


ization ( Iliev et al.| 20141 in a box of comoving length 
425/i“^ Mpc (corresponding to ~ 4 deg at z=7) with 
165 billion particles distributed on a grid of 10, 976^ cells 
(3.9 h~^kpc gravity force resolution) and a radiative transfer 
grid of 504^ cells. The N-body simulation has been run from 
redshift z = 300 to 2: = 2.6, with initial conditions generated 
using the Zel’dovich approximation and a power spectrum 


of the linear fluctuations given by the GAME code ( Lewis, 
Challinor fc Lasenby|[2000 1 . This simulation was then used 
as an input to the radiative transfer code C^-RAY (Mellema 
et al.|[^06a l to follow the reionization history of the IGM. 
More specifically, the halo catalogues were used to construct 


the sources of ionizing radiation as in Iliev et al. (20071. 
As the minimum resolved halo mass is Mh,min = 10“ Mq, 
haloes with masses of 10® — 10 ® Mq were modeled as a sub¬ 


grid population (Iliev et al. 2014 Ahn et al. 20151. All haloes 
were assigned an ionizing photon emission rate per unit time, 
Nry, proportional to the halo mass, M^: 


^ _ gjMh0. b f 10 Myr 


At 


( 1 ) 


where nip is the proton mass, D(, and Do have their usual 
cosmological meaning. At = 11.46 Myr is the time between 
two snapshots of the N-body simulation, and g-y is a source 
efficiency coefficient that incorporates the star formation ef¬ 
ficiency, the total photon production per stellar baryon and 
the io nizing photon escape fraction (Iliev et al.|2006 Jensen 


et al. 


20141. Haloes with masses down to 10“ Mq were 


assigned a source efficiency of g-y = 1.7. Smaller sources 
with masses down to 10® Mq were assigned g^ = 7.1, to 
account for a lower metallicity and a more top-heavy ini¬ 
tial mass function, but they were assumed to be suppressed 
within ionized regions (for ionization fraction higher than 
10 per cent; Iliev et al.|2014 (. The radiation emitted by the 
sources is propagated through the gridded density field, and 
the distribution of neutral hydrogen is obtained at various 
redshifts. This is used to calculate the associated differen¬ 
tial brightness temperature according to the usual formalism 
(e.g. Field|[l959 Madau, Meiksin fc Rees||1997 Furlanetto, 


Oh fc Briggs|2006D : 


STb = 28.5 mK (1 -f S)xhi ( 

\ ^ j HI 1 ^ 0.042 0 . 73 / 


l-hz 

10 


0.24 \ 
D^ j 


1/2 


( 2 ) 


where = nni/{nn) is the mean density of neutral 

hydrogen in units of the mean density of hydrogen at redshift 

2. 

For our purposes, we used boxes from the simulation 
at 2 =6.68, 7.06 and 7.3, corresponding to volume (mass) 
averaged ionized fractions {x) = 0.93 (0.95), 0.65 (0.73) and 
0.48 (0.58), respectively. These particular boxes were chosen 
because HSC will have two narrow-band filters observing at 
redshifts 6.6 and 7.3, while 7.06 is an intermediate value. 

The same simulations were processed with a Lya ra¬ 
diative transfer code to model high-2 LAEs and study their 
observability. Motivated by detailed radiative transfer cal¬ 
culations by |Laursen, Sommer-Larsen fc Razoumo^ ( |2011[ ), 
the Lya line was modeled as a double peaked profile with lit¬ 
tle emission at the line centre, and a width that depends on 
the halo mass. Intrinsic luminosities were calibrated against 
observations, with a model where the Lya luminosities of 
haloes of a given mass follow a log-normal distribution with 
a mean that is proportional to the halo mass. After assigning 
an intrinsic Lya spectrum to the dark matter haloes in the 
N-body simulations, the observed luminosities are calculated 
including the attenuation from the IGM along a large num¬ 
ber of lines of sight from each of the haloes (for more details 
we refer the reader to the original papers [jensen et al.|2013'l 
20141. From the same work we extracted the Lya intrinsic 


and transmitted luminosities, which we use to produce HSC 
mock observations. 

In the computation of the 3D cross-power spectra we 
merged bins to obtain AA: > 0.02 h Mpc which corre¬ 
sponds to the smallest mode resolved by a FoV of 16 deg®, 
i.e. equivalent to our simulations. We also made sure to avoid 
correlations in power due to the window functiorQby using 
a binning with A log k = 0.02. As the FoV used to com¬ 
pute the 2D cross-power spectra is smaller (i.e. 7 deg® and 
1.7 deg® at 2 =6.6 and 7.3, respectively; see section 5), we 
used Alogfe = 0.03 (0.05) and AA: > 0.04 (0.07) h Mpc 
for 2 = 6.6 (7.3). 


3 CROSS-POWER SPECTRUM 

In this section we present our calculations of the theoretical 
and observational cross-power spectra. 

At each redshift, the 21cm-galaxy cross-power spectrum 
at wave number A: = |A:|, A^i^ga/A:), can be decomposed into 


^ The sphere used to compute a spherically averaged P{k) in a 
simulation of comoving length 425 h“^Mpc, must be equivalent 
in volume and thus have a radius R = 264 /i“^Mpc comoving. A 
window function for a spherical tophat has its first zero at dfc ■ ij ~ 
4.5, so that fc-values spaced by less than 4.5/A? = 0.02 h Mpc“^ 
will be correlated ([Feldman, Kaiser fc Peacock||1994| | Furlanetto] 


&: Lidz|2007 Lidz et al.|2009 Wiersma et al.|2013i. 



























































4 D. Vrhanec et al. 


three contributing terms (e.g. |Lidz et al.|2009 l: 

^21,gal{k) = A2l^gal{k)/STl,0 

= {xhi)[A^^j g^i{k) + Apg^i(k) + A^^jp g^i{k)], (3) 

where Al^^^g^i, and are the neutral 

fraction-galaxy, density-galaxy and neutral density-galaxy 
cross-power spectra, respectively. STi,q is the 21cm bright¬ 
ness temperature relative to the CMB for a fully neutral 
gas element at the mean cosmic density, and (xhi) is the 
volume-averaged neutral hydrogen fraction. A^ is the di¬ 
mensionless cross-power spectrum of two random helds, a 
and 6, and it is equal to A^ j,(fc) = k?'Pa,b{k)/2 tx^ for the 3D 
cross-power spectrum, and A^ i,(fc) = 2'Kk^Pa,b{k) for the 
2D power spectrum. Pa,6 represents the dimensional cross¬ 
power spectrum between helds a and b. The latter are repre¬ 
sented in terms of their fractional huctuations at a location 
r, i.e. Sa{r) = {a{r) — (a))/(a), and similarly for lP| A more 
detailed discussion of the various terms can be found in lLidzl 

[eTaTIpOO^ . 

3.1 Theoretical 21cin-galaxy cross-power 
spectrum 

To understand the 21cm-galaxy cross-power spectrum we 
hrst show the theoretical spherically averaged 3D 21cm- 
galaxy cross-power spectrum at z = 7.3 (which was com¬ 
puted using haloes with Mji > 10^° Mq, i.e. 3 million 
galaxies at z = 6.68, 2.3 million at 2 : = 7.06 and 1.9 
million at 2 = 7.3), together with its contributing terms 
(Fig-iii top panel) and the corresponding cross-correlation 
coefficient (Fig. bottom panel), defined as r 2 i,gai{k) = 
P 2 i,gai{k)/[P 2 i{k)Pgai{k)Y^'^ . This correspouds to the ideal 
case in which all galaxies could be observed. From the be¬ 
haviour of A^p gai it is clear that, as expected, the galax¬ 
ies are strongly correlated with the density held on small 
scales, because galaxy formation is biased toward high den¬ 
sity regions, while the correlation decreases as we move to¬ 
wards larger scales, but always remains positive. The neutral 
hydrogen-galaxy cross-power spectrum, A^^^ gai, instead, is 
negative on large scales where there is a paucity of galaxies 
but most of the HI resides. A turn around is observed in cor¬ 
respondence of the typical scale of the HII regions, and then 
the correlation drops off since the hydrogen inside such re¬ 
gions is completely ionized independently from the number 
of sources. A^^^p gai is positive on the largest scales and be¬ 
comes negative towards smaller scales, where it cancels out 
with A^p gai- The hnal 21cm-galaxy cross-power spectrum 
thus follows the shape of A^^^ gai on small scales, and that 
of A^^^ gai and A^ gai on large scales. The cross-correlation 
coefficient (bottom panel of Fig. shows more clearly that 
the 21cm signal and the high -2 galaxies are anti-correlated 


® Note that we evaluate the theoretical cross-power spectrum 
with (a) = o,i)/A, where N is the number of pixels in the 

portion of the simulation used. All the quantities are calculated 
like this, with the exception of the galaxy field in mock observa¬ 
tions, which is instead calculated using {Ngai) = Ngai/V, where 
Agal is the number of galaxies in the mock observation, and V is 
the volume of the survey. This was done for an easier comparison 
with the shot noise power spectrum Pghot{k) = l/ugaj, where 
ngai is the average number of galaxies in the survey volume. 



Figure 1. Top panel: spherically averaged 3D 21cm-galaxy 
cross-power spectrum, (black solid line) at 2 = 7.3, to¬ 

gether with its contributing terms, i.e. the neutral hydrogen- 
galaxy cross-power spectrum, (blue dashed), the 

density-galaxy cross-power spectrum, (green dashed- 

dotted) and the neutral density-galaxy cross-power spectrum, 
p gal dotted). Bottom panel: 21cm-galaxy cross¬ 

correlation coefficient, r 2 i,gal (black solid), and zero-correlation 
coefficient (black dotted). 


on large scales, and become uncorrelated on scales smaller 
then the typical size of the ionized regions. 

Similar conclusions were drawn by Lidz et al. (20091 


and Wiersma et al. (20131, although our results are closer 
to those of|Lidz et ^ (20091 because of the lower resolution 


employed in the simulation by Wiersma et al. (20131. 

Figure shows A^i^gai s-tid ^ 21 ,gai for the chosen red- 
shifts. We can see that the amplitude of the power spectrum 
decreases with decreasing redshift, while the turnover point 
shifts towards larger scales. This indicates that, as reioniza¬ 
tion proceeds, the anti-correlation decreases because of the 
paucity of neutral hydrogen, and the ionized bubbles grow 
in size. This is more clearly seen in the behaviour of the 
cross-correlation coefficients, which shift towards smaller k 
with decreasing redshift. 


3.2 Observed 21cm-LAE cross-power spectrum 

In this section we will show our predictions for the 2D 21cm- 
LAE cross-power spectrum as it would be observed with 
LOFAR and HSC. To do that, we added projection effects 
and constrained the galaxy number density to match HSC 
expectations, and we added noise to the 21cm field to sim¬ 
ulate LOFAR observations. 

HSC will probe the reionization epoch with the Ultra¬ 
deep and Deep layers of the HSC Survey. Observations 
are made with narrow-band filters (A 2 = 0.1, equiva¬ 
lent to approximately one tenth of our simulation length 
~ 42 /i~^Mpc), so that the LAEs redshift will be tightly 
constrained. Because the LAEs detected with a particular 
filter will be observed as if they were lying on a single plane, 
the observed 21cm-LAE cross-power spectrum will be a cir¬ 
cularly averaged 2D cross-power spectrum. HSC will observe 
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Figure 2. Top panel: spherically averaged 3D 21cm-galaxy 
cross-power spectrum, at 2 = 7.3 (black solid line), 7.06 

(green dashed), and 6.68 (red dashed-dotted). Bottom panel: 
21cm-galaxy cross-correlation coefficient, r 2 i,gai-, corresponding 



Figure 3. Top panels: 2D unnormalized, circularly averaged 
21cm-galaxy cross-power spectra at 2 : = 7.3 (left panel) and 
6.6 (right). Shaded areas indicate scatter from 10 mock observa¬ 
tions. Bottom panels: 21cm-galaxy cross-correlation coefficient, 
r 2 l,gal, corresponding to 


4 fields of 7 deg^ at redshift z = 6.6 as part of a Deep layer, 
and 4 fields of 1.7 deg^ (two at z = 6.6 and two at z = 7.3) 
as part of a Ultra-deep layer (Ouchi 2012, private commu¬ 
nication). One of the fields in the Deep layer is ELAIS-Nl, 
which will also be observed with LOFAR ( Jelic et al.|2014 1 . 

We reduced the box dimension to match the HSC’s field 
of view (7 deg^ at ^ = 6.6 and 1.7 deg^ at 2 = 7.3) by remov¬ 
ing external cellf^ We then divided our simulation boxes of 
brightness temperature and galaxies into 10 sub-boxes of 50 
slices each, corresponding to a A2 = 0.1. Each sub-box ob¬ 
tained from the galaxy simulation is collapsed onto a single 
plane to mimic the fact that HSC observations will provide 
a 2D map of galaxies. This map is then correlated with each 
of the 50 slices of the corresponding brightness temperature 
sub-box to obtain 50 2D 21cm-galaxy cross-power spectra, 
which are then averaged to mimic the result of observations 
of a single FoV. From the 10 sub-boxes we then obtain 10 
2D 21cm-galaxy cross-power spectra, which can again be 
averaged so that our results are not sample dependent. 

Figure 1^ shows final, unnormalized by STbo, 2D 21cm- 
galaxy cross-power spectra before including the noise and 
the constraints on the galaxy number densitj*^ Even in 
2D the cross-power spectra still retain much of their shape, 
although some features are lost due to projection effects and 
reduction in field of view, e.g. the turnover point is not clear 
anymore. Projection effects also induce a reduction in the 
value of the anti-correlation, clearly observed in the cross- 


® The choice of removing external cells is arbitrary and we have 
checked that it does not affect the final results. 

We note that while the solid lines represent the absolute value 
of the average cross-power spectrum (i.e. the average could be 
both positive and negative), the shaded area is obtained from the 
scatter in absolute averaged values (i.e. only positive numbers). 
For this reason the solid lines do not always lie at the center of 
the shaded areas. 


correlation coefficient, which drops from r 2 i,gai « — 1 to 
r 21 ,gal ~ -0.5. 

When selecting LAEs for our mock observations we as¬ 
sign intrinsic equivalent widths (EW) to the LAE sample ac¬ 
cording to a log-normal distribution, as was done by [Jensen] 
et al. (20141. The distribution is designed to approximately 
fit observations made by Jiang et al. (20131, while giving 
65 % of the galaxies EW below 25 A (consistent withjStark 


et al. 2010] as shown in Fig. 1 in Jensen et al. 20141. We 


first selected all the galaxies with EW > 20 A (consistent 
with HSC expectations), and among these only the 1375 (20) 
most luminous ones at 2 = 6.6 (7.3), to match the number 
expected to be observed by HSC. 

LOFAR will be detecting the cosmological 21cm signal 
with a field of view of 5 x 5 deg^, and an angular resolu¬ 
tion of 3.5' ( Zaroubi et al.|[2012 |. To simulate the LOFAR 
noise at each frequency, we filled a LOFAR measurement 
set (the real and imaginary parts of the visibilities) with 
Gaussian random numbers. This was then imaged (Fourier 
transformed, accounting for the proper weighting) to obtain 
noise maps in real-space, and their root mean square was 
normalised according to (e.g. jTaylor, Carilli fc Per ley] 1999] ): 

W SEFD 

— / , ' ( 4 ) 

Vs ^2N{N - l)AvtM 

where W is a factor used to increase the noise according to 
the adopted weighting scheme, Vs is the system efficiency, 
SEFD is the system equivalent flux density, N is the num¬ 
ber of stations, Au is the bandwidth, and tint is the inte¬ 
gration time. Based on empirical SEFD values for LOFAR 
(e.g. 3000 Jy at 150 MHz towards the zenith; jvan Haai^ 
lem et al. 20131, we expect (t„ to be about 76 mK at a 
resolution of 3.5 arcmin, at 150 MHz, after 600 hours and 
0.5 MHz of integration and assuming N = 48, W = 1.3, 
Vs = 0.9. Note that adopted noise values are indicative only, 
and they may change in the actual observations due to e.g. 
time-variable station projection losses of sensitivity, smaller 


system efficiency, etc. (van Haarlem et al. 20131. More de- 
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Figure 4. Top panels: 2D unnormalized, circularly averaged 
21cm-LAE cross-power spectra at 2 = 7.3 (left panel) and 6.6 
(right). Shaded areas indicate scatter from 10 mock observa¬ 
tions. Solid (dashed) lines refer to the cross-power spectra with 
(without) LOFAR noise. Bottom panels: 21cm-galaxy cross¬ 
correlation coefficient, r 2 i,gai, corresponding to 


tails about simulating the LOFAR noise can be found in e.g. 
Patil et al.| ( |2014| ). The simulated LOFAR noise was added 
to the brightness temperature map from the simulation. 

In Figure]^ we plot the resulting 2D unnormalized, cir¬ 
cularly averaged 21cm-LAE cross-power spectra with (solid 
lines) and without (dashed) LOFAR noise. Despite the spec¬ 
tra being much noisier than the previous ones at all scales, 
a dependence of the normalization on redshift (i.e. amount 
of HI) and an anti-correlation {r 2 i,gai ~ —0.20) are still vis¬ 
ible on large scales, with levels of slgnihcance of p = 0.04 
at z = 6.6 and p = 0.048 at 2 = 7.3, although the turnover 
point can not be clearly identified. Even without LOFAR 
noise, observations at small scales will still be largely af¬ 
fected by shot noise and will not offer any reliable data (see 
sec.[^. From this analysis we conclude that only scales larger 
than ~ 60 (45) h~^ Mpc, i.e. k < 0.1 (0.14) h Mpc“^, at 
2 =6.6 (7.3) can be used for cross-correlation studies. 

For a HSC FoV equal to the one of LOFAR, though, we 
would expect the detection of 3140 and 90 LAEs at 2 = 6.6 
and 7.3, respectively. In this case (see Eig. the overall 
noise would be reduced, the anti-correlation signal would 
be stronger {r 2 \,gai ~ —0.30), large scales could be more 
reliably used, and information could be extracted down to 
~ 60 (30) h~^ Mpc, i.e. k > 0.1 (0.2) h Mpc“^, at 2 = 6.6 
(7.3). In addition, also information at scales larger than ~ 
130 (80) h~^ Mpc at 2 = 6.6 (7.3) and up to ~ 310 h~^ Mpc 
would be available. 


4 CROSS-CORRELATION FUNCTION 

In this section we present our calculations of the theoretical 
and observational cross-correlation functions. 



Figure 5. Same as Figure [4] but for a HSC field of view of 16 
deg^ at both redshifts. 

4.1 Theoretical 21cm-galaxy cross-correlation 
function 

In addition to the cross-power spectrum, we have also com¬ 
puted the cross-correlation function, which shows how the 
correlation between two fields changes in real space. The 
cross-correlation function between fields a and b is defined 
as ^a,i)(r) = (5a(x)(5i,(x -|- r)), where (5(x) is the fractional 
fluctuation of the held at location x. 

The 3D 21cm-galaxy cross-correlation function can then 
be calculated from the cross-power spectrum as (|Park et al.| 
2014p ] 

^21.gal{r) = I P21,gal{kf-^A-Kk^dk. (5) 

In Figure we show the 21cm-galaxy cross-correlation 
function at 2 = 7.3, together with the different terms that 
contribute to it. ^p,gai shows positive correlation on small 
scales, where there is an overdensity of both gas and galaxies, 
and no correlation on large scales. The neutral and galaxy 
helds, ^XHi,gai, are anti-correlated on small scales (where the 
gas is mostly ionized and there is an overdensity of galaxies), 
mildly correlated on scales just larger than the typical scale 
of ionized bubbles (where neutral hydrogen is more abun¬ 
dant), and show no correlation on large scales (where most 
of the neutral hydrogen resides, but there is a paucity of 
galaxies). E,XHiP,go.i and (,21,gai behave similarly to ^^Hj.gai, 
although ^XHgp,gai turns over to positive values and no cor¬ 
relation on much smaller scales. We can see that the typical 
scale of ionized regions is ~ 50 h~^ Mpc. 

In Figure we show the theoretical 3D 21cm-galaxy 
cross-correlation function at 2 = 6.68, 7.06 and 7.3. The 

Note that when computing the cross-correlation function from 
the cross-power spectrum, uncertainties arise because of the inte¬ 
gration over a finite box size and finite resolution (e.g. uncertain¬ 
ties in the information about the turn over scale; [^ark et al.|2014[ |. 
However, because of the large box size and number of galaxies, it 
is computationally much more efficient to compute the 3D 21cm- 
galaxy cross-correlation function from the cross-power spectrum 
than directly. 
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Figure 6. Theoretical 3D cross-correlation functions at z = 
7.3 multiplied by 10 for better resolution: 21cm-galaxy, ^21,gal 
(black solid line), neutral fraction-galaxy, ixn 1 , 90.1 (blue dashed), 
density-galaxy, ^p,gai (green dashed-dotted) and neutral density- 
galaxy, ixHip, 9 oi dotted). The black dotted line indicates 

zero correlation. 



r (h ^ Mpc) r (h ^ Mpc) 


Figure 8. 2D 21cm-LAE cross-correlation function for our refer¬ 
ence mock observations at 2 = 7.3 (top left panel) and 6.6 (top 
right) and for mock observations with fields of view of 16 deg^ 
at 2 = 7.3 (bottom left) and 2 = 6.6 (bottom right). The black 
dotted lines indicate zero correlation and shaded areas indicate 
scatter from 10 mock observations. 



Figure 7. Theoretical 3D 21cm-galaxy cross-correlation function 
at 2 = 6.68 (red dashed-dotted line), 7.06 (green dashed) and 7.3 
(black). The black dotted line indicates zero correlation. 


qualitative behaviour of the curves is similar, with an anti¬ 
correlation on small scales, indicating the typical scale of the 
ionized regions, followed by a small positive correlation, and 
no correlation on larger scales. As for the case of the power 
spectrum, the anti-correlation is smaller with decreasing red- 
shift due to the fainter 21cm signal. 


4.2 Observed 21cm-LAE cross-correlation 
function 


The observed 2D cross-correlation function can be calcu¬ 
lated as: 


^21,gal{r) 


Slae{x)52i{x + r) 
^pair {r ) 


( 6 ) 


where and ^21 are fractional fluctuations of the LAE 

and 21cni fields, respectively, and Npair(r) is the number of 
21cm-LAE pairs at a separation r. 

In Figurej^we plot the 2D 21cm-LAE cross-correlation 
functions for our reference mock observations (i.e. the equiv¬ 
alent of Fig. 1^, together with those expected for a larger 
held of view of 16 deg^ (i.e. the equivalent of Fig. |^. The 
observed cross-correlation functions show a behavior similar 
to the theoretical ones. Noise is large at all scales, resulting 
in a large scatter. While the average of 10 mock observations 
for both redshifts shows clear anti-correlation at small scales 
which goes towards no correlation at large scales, scatter is 
large, so the detection of the anti-correlation might not be 
possible in a single mock observation. The anti-correlations 
become much clearer in larger Helds of view, especially at 
redshift 6.6. 


5 DISCUSSION 

Observations of 21cm emission and high-a LAEs are ex¬ 
tremely challenging, and both will suffer from severe noise 
problems. Even assuming that foregrounds subtraction will 
work perfectly, the system noise will still largely exceed the 
expected signal, in particular at the smaller scales, so that 
possibly only scales larger than ~ 60 h~^ Mpc (correspond¬ 
ing to fc ~ O.lh Mpc“^) will be accessible by a telescope 
like LOFAR. In addition, the Held of view of HSC is much 
smaller than that of LOFAR, so that only a fraction of the 
large scales observed by LOFAR will be covered also by HSC. 

To illustrate this issue further, in Figure]^ we show the 
2D 21cm auto-power spectra with and without the LOFAR 
noise after 600 hours of observation in Helds of view of 7 
deg^ at 2 = 6.6 and 1.7 deg^ at 2 = 7.3, i.e. equivalent to 
the ones of hscE!1 At both redshifts noise on scales smaller 

Note that in observations of the 21cm auto-power spectra the 
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Figure 9. 2D 21cm auto-power spectrum with (solid lines) and Figure 10. 2D LAE auto-power spectra (solid lines) and shot 

without (dashed) LOFAR noise at z = 6.6 (red) and 7.3 (black). noise power spectra (dashed) at z = 6.6 (red) and 7.3 (black). 


than ~ 60h~^ Mpc is orders of magnitude larger than the 
expected signal, while it decreases gradually at larger scales. 
Noise on large scales at z = 6.6 is somewhat larger than at 
2 = 7.3. 

The HSC observations discussed here are groundbreak¬ 
ing, as they will increase the number of detected high-z 
LAEs by at least one order of magnitude. However, substan¬ 
tial shot noise is still expected, as shown in Figure Ob¬ 
servations at both redshifts will be dominated by shot noise 
at scales below 10h~^ Mpc (i.e. k > 0.6h Mpc“^) at z=6.6, 
and Mpc (i.e. k > 0.2h Mpc“^) at z=7.3, while at 

large scales the LAEs auto-power spectrum is stronger than 
that of the shot noise, in particular at z = 6.6. 

Since the 21cm-LAE cross-power spectrum will be af¬ 
fected by noise from both instruments, we expect to be 
able to probe only scales larger than ~ 60/i“^ Mpc (i.e. 
k < O.lh Mpc“^). Such scales will still have shot noise, in 
particular at z = 7.3, but this should not prevent the detec¬ 
tion of an anti-correlation. 

Stronger anti-correlation could be detected by reduc¬ 
ing the noise, e.g. increasing the integration time for 21cm 
observations {anoise ~ LOFAR; Eq. |^, or with a 

larger field of view. The latter would Increase the number of 
observed LAEs and thus reduce the shot noise and extend 
the number of observed fc-modes. 

In 21cm-LAE cross-correlation function noise is large 
at all scales. This is because, unlike in cross-power spec¬ 
trum, noise does not get separated by its fc-modes, and 
thus it is equally distributed on all scales. Comparing cross¬ 
correlation functions at z = 7.3 for our reference mock ob¬ 
servations and for mock observations with a larger held of 
view (Fig. |^, in the latter case a larger amplitude of the 
anti-correlation as well as a smaller scatter can be observed 
because of a reduction of the noise component. While the 


expectation value of the noise power spectrum can be subtracted 
from the measurements. However, in observations of the 21cm- 
LAE cross-power spectra this is not possible, since instrumental 
effects from LAE observations are also present and their influence 
can not be treated separately. The same reasoning applies to LAE 
observations. 


scatter is smaller also at z = 6.6, the amplitude is not in¬ 
creased. We suggest that this is due to the LOFAR noise, 
which has a stronger effect at z = 6.6, despite being smaller 
in absolute terms. The noise should also be responsible for 
the positive correlation observed at large scales in the top 
panels of Figures EH and |12| 

To reduce the noise levels we also smoothed the 21cm 
held with a Gaussian hlter with standard deviations of 
cr = 2, 5 and 10 (with smoothing radii 4.14 h~^ Mpc, 
10.35 h~^ Mpc and 20.70 h~^ Mpc, respectively) at z = 6.6 
(Fig. 111 and ct = 1, 2 and 5 (with smoothing radii 2.07/i“^ 
Mpc, 4.14 h~^ Mpc and 10.35 h~^ Mpc, respectively) at 
z = 7.3 (Fig. |12| l. We clearly see that both the average 
cross-correlation function and the scatter become smoother 
with increasing a. However, by smoothing the held we also 
loose information (e.g. in terms of anti-correlation ampli¬ 
tude), which is visible when comparing results with different 
a. Smoothing the signal at z = 6.6 reduces the noise on small 
scales enough that the anti-correlation becomes clear even 
for the largest scatter. At z = 7.3, instead, the shot noise 
is larger because of the smaller LAE sample, so that even 
after smoothing the scatter on small scales remains large. 

While smoothing reduces noise in the cross-correlation 
function (which remains though still noisy on all scales), it 
is not helpful when applied to the 21cm-galaxy cross-power 
spectrum. The reason for this is that because the noise is be¬ 
ing separated by its fc-modes, smoothing would affect only 
the small scales which would still be over-contaminated by 
shot noise. However, the separation of noise by its fc-modes 
is exactly what makes large scales observable and the 21cm- 
galaxy cross-power spectrum a more useful probe of reion- 
Ization. 


6 CONCLUSIONS 


In this paper we have investigated the properties of the 
21cm-galaxy cross-power spectrum and its observability 
with upcoming observations by LOFAR and HSC. To this 
aim, we have used snapshots at redshifts 6.68, 7.06 and 7.3 
from N-body + radiative transfer simulations (Iliev et al. 
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Figure 11. 2D 21cm-LAE cross-correlation function at 2 = 6.6 
(top left panel) using a Gaussian filter with standard deviation 
(T = 2 (top right), cr = 5 (bottom left) and ct = 10 (bottom right). 
The black dotted lines indicate zero correlation and shaded areas 
scatter from 10 mock observations. 
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Figure 12. 2D 21cm-LAE cross-correlation function at 2 = 7.3 
(top left panel) using a Gaussian filter with standard deviation 
(T = 1 (top right), fT = 2 (bottom left) and a = 5 (bottom right). 
The black dotted lines indicate zero correlation and shaded areas 
scatter from 10 mock observations. 


2014), processed with a Lya radiative transfer code (Jensen 


et al. 20131. Our theoretical 3D 21cm-galaxy cross-power 
spectrum agrees with previous investigations (i.e.|Lidz et al.| 


2009 Wiersma et al. 2013). More specifically, we are able 


to recover the same redshift dependence and shape, with a 
distinct turnover point indicating the typical scale of ionized 
bubbles. We confirm that the 21cm-galaxy cross-power spec¬ 
trum could provide information on the progress of reioniza¬ 
tion and the typical size of HIT regions at different redshifts. 

The measured 21cm-LAE cross-power spectrum suffers 
from projection effects (as it is 2D), as well as from noise in 
both radio and LAEs detections. LOFAR recently started 
observations of 21cm emission from neutral hydrogen in the 


redshift range 2 = 6—11.4 (Yatawatta et al. 2013 Jelic 


et al. 2014), while HSC will also soon start its observa¬ 


tional campaign with two narrow-band filters searching for 
LAEs at 2 = 6.6 and 7.3 (M. Ouchi 2012, private conversa¬ 
tion). Both telescopes plan to observe the ELAIS-Nl field 
at 2 = 6.6, making it possible to detect the 21cm-galaxy 
cross-power spectrum. We constructed mock observations 
specifically tailored to match LOFAR and HSC campaigns 
at redshifts 6.6 and 7.3. Our mock observations show that 
despite the observed spectra being much noisier than the 
corresponding theoretical 3D ones, dependence of the nor¬ 
malization on redshift (i.e. amount of HI) is clearly visible, 
as well as the anti-correlation between the two fields, with a 
cross-correlation coefficient r 2 \,gai ~ —0.20 at levels of sig¬ 
nificance of p = 0.04 at 2 = 6.6 and p = 0.048 at 2 = 7.3. 
However, the turnover point can not be clearly determined 
because small scales will be overwhelmed by noise. 

We also investigated properties and observability of the 
21 cm-galaxy cross-correlation functions, which are expected 
to be negative on small scales, mildly correlated on scales 
just larger than the typical size of ionized regions, and show 
no correlation on even larger scales. This agrees well with 
predictions of the 21cm-galaxy cross-correlation function by 


Park et al. (2014). Despite observational effects like noise 


and galaxy number densities, observed correlation functions 
should retain the theoretical shape. However, unlike the 
observed 21cm-LAE cross-power spectrum, the correlation 
function suffers from a strong noise on all scales (as it does 
not get separated by its fc-modes), thus uncertainties in the 
signal will be large. 

In summary, the 21cm-LAE cross-power spectrum is a 
powerful probe of the EoR which could provide invaluable 
information on the progress of reionization and the typical 
scale of ionized regions at different redshifts. Observations 
with LOFAR and HSC will finally make detection of the 
21cm-LAE cross-power spectrum possible at redshift 6.6, as 
they both plan to observe the ELAIS-Nl field. These obser¬ 
vations are going to be very challenging and have substantial 
problems with noise, but they will still be able to detect the 
large scales of the cross-power spectrum, which is expected 
to show an anti-correlation between the two fields. 
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